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ABSTRACT 

A derailed gravimetric geoid for the United States, 
which has been computed from a combination of 
satellite derived gravity data and local gravimetric 
measurements, is presented herein. The detailed 
geoid is shown to possess a precision of ±2 meters 
rms. Comparisons of the detailed geoid with results 
of Rapp, Fischer, and Rice are presented. 
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A DETAILED GRAVIMETRIC GEOID FOR 
THE UNITED STATES 


1. INTRODUCTION 

The objective of the work described within this report was to compute a 
gravimetric geoid for the United States having an rms accuracy of ±2 meters. 
This objective was achieved through the use of a combination of surface and 
satellite gravimetric information. The detailed geoid presented here is believed 
to represent the first detailed gravimetric geoid of the entire United States ever 
published in the open literature. Rapp (1967) has previously published a com- 
bination geoid for a portion of the United States. 

Comparisons of the computed gravimetric geoid with the results of Rapp, 
and with the astrogeodetic geoids of Fischer (1967) and Rice (1962, 1970), indi- 
cate that the anticipated ±2 meter rms accuracy has been achieved. 

2. METHOD OF COMPUTATION 

The geoidal undulation at any point P on the earth can be computed using 
the well known Stokes' formula: 

2 it ir/2 

N(c P’ X)= 4^T [ f Ag T (cp , ,X / ) S('I') cos cp'dcp'dx/ (1) 

/ 1 v 

X =0 cp = - 2 


1 


where: 


<P , K = The latitude and longitude, respectively, of the computation point. 

<P ‘ . X ' = The latitude and longitude, respectively, of the variable integration 
point. 

N (<p , \) - Geoid undulation at <P , X . 

R = Mean radius of the earth. 

G = Product of the universal gravitational constant and the mass of the 
earth. 

Ag T (<p',\') c Free air gravity anomaly at the variable points' , X' . 

S(V) = 6 sin (V/2) + 1 - 5 cos Y 

V ' sin(4»/2) V 

- 3 cos 4* In (sin (V/2) + sin 2 (V/2) 

where 


4* = cos -1 [sin cp sin cp’+ cos<p cos5p' cos (\ - X*)] 

In order to combine surface and satellite gravity data for geoid computation 
the earth is divided into two areas, a local area (A j ) surrounding the point P, 
and the remainder of the earth (A 2 ). Also the anomalous gravity in each area 
is partitioned into two parts represented by the symbols Ag s and Ag 2 . The Ag s 
values are defined as that part of the anomalous gravity field which can be 
represented by the coefficients in a satellite derived spherical harmonic expan- 
sion of the gravitation potential. The 1969 SAO Standard Earth given in Table 1 




was used in all computations described in this paper. The Ag 2 values are de- 
fined as the remainder of the anomalous gravity field. Using this division of the 
earths surface in two areas and of the anomalous gravity each into two components 
one can write equation (1) in the form: 


where 


N 


l 


R 

4ttG 



N(cp, X) = Nj + N 2 + N 3 


(2) 



Ag g (cp', X ') S (ty) cos cp' dcp' d\ ' 


Ag 2 (cp ' » X ' ) S (ty) cos cp' dcp ' dX ' 


( 3 ) 




R 

4ffG 



ag 2 (cp', X ') S ('I') cos cp' dcp' dX ' 


The following paragraphs discuss how each of the three components presented 
in equation (3) is handled in the computations. 


Given a set of satellite derived coefficients in the spherical harmonic ex- 
pansion of the gravitational potential a number of methods exist for computation 
of the N t component of the geoid undulation. 
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Table 1 

Normalized Coefficients for the SAO-69 Standard Earth Model 


(Gaposchkin and Lambeck, 1970) 
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The computation of Nj was not carried out in the present case by using the 
integration indicated in equation (3). Rather the procedure described by Bacon, 
et al., (1970), was used. Briefly this procedure consists of fixing a value of the 
potential, W Q , and computing the component N 4 as 


=R g -R e 


where: 


R G is the radial distance of the equipotential surface defined by W o and 
the potential coefficients of the SAO 1969 Standard Earth. 

R is the radial distance to a selected reference ellipsoid defined by a 
semimajor axis (a ) and flattening (f). 


The radial distance, R G , to the equipotential surface W o at a particular 
latitude and longitude , \ t is determined by using the equation 


k m . , n 

W = } } (~ J {C sin m\. + S sin (5) 

o / t / t \ f J nm * nm 1 ' ' 


n = 2 m a 1 


P (<M 

nm v * ' 


The only unknown in this equation is r. Using an iterative three point in- 
verse interpolation scheme the value of r (i.e., r = R ) which will make equation 

G 

5 an identity is determined. Using this value of R c , and R £ computed using 
the input values of a and f of the reference ellipsoid, a geoid undulation com- 
ponent N , is computed. 
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For the computations described in this paper the area A , for a point at 
which the geoid was being computed was defined to consist of a twenty degree 
by twenty degree area centered on the computation point. The computational 
formula used was: 

400 

N 2 = — ... ^ A «2 <vi> kj) S OP,) cos <pj A 9P* d ( 6 ) 

n U J-i 

where 

A g 2 (qpf , \* ) is the mean value of A g 2 within the J th 1° x 1° square 
S(*P, ) is the value of Stokes* function at the center of the j th 1° x 1° square. 

Af' = A \'= 1°. 

The value of Ag 2 used for each 1° x 1° square was computed using the 
formula 


A g 2 = A g e - A g f 

The ~Kg e values are mean 1° x 1° free-air anomalies provided by surface 
gravity data. The primary source material for this gravity data is indicated in 
Appendix A. 

The Ag s values are that part of the mean 1° x 1° free-air anomalies repre- 
sented by the satellite harmonic coefficients used in computing N, . The A g, 
values are computed using the formula 


6 




k n 

^e=>e/ / ( n " 1) CO,,,, cos m + S sinmV] P (<p"\ 

c r / / , nm nm nm V ^ ' 


n*2 m a 0 


where 

7 e = Equatorial gravity in mllligals (0.98 x 10 +6 ). 

C nm , S nm = Normalized geopotential spherical harmonic coefficients 
C ,S except for C and C . 
k = Upper limit on degree and order of the geopotential model, 
n = Degree index of harmonic coefficients, 
m - Order index of harmonic coefficients. 

In equation (5), the C 2Q and C 4Q terms do not represent the complete co- 
efficients but rather the difference between the complete coefficients and the 
coefficients compatible with the ellipsoid used in computing N } . In order for 
the above described procedure to produce correct results, the quantities A"g e , 

A g g , and the a and f which define the ellipsoid used to compute N j must all be 
compatible. Compatibility implies that the values of C 2Q and C 4Q used to com- 
pute the values of theoretical gravity needed to obtain Ag e and Ag, are the same 
as the values of C 2Q and C 4Q implied by the reference ellipsoid. Correct re- 
sults in absolute sense are also dependent upon the value of W Q being chosen 
to represent the true value of the potential of the geoid. The effects of not mak- 
ing a, and f compatible are twofold. First, all the computed geoid 

heights may by in error by a constant; in addition, there will be a systematic 
error as a function of latitude. The effect of selecting an incorrect value of W q 
would be to introduce a constant error in all geoid heights. 
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In the calculations described here the term N 3 in equation ( 2 ) is set equal 
to zero. This is equivalent to assuming that the satellite derived approximation 
to the gravity field is adequate for the area A , at a distance of greater than ten 
degreees from the computation point. 

3. FINAL RESULTS 

t 

The final detailed gravimetric geoid is presented in Figure 1 using a one 
meter contour interval. The geoid is referenced to an ellipsoid with a flattening 
given by f = 1/298.255. 

The primary sources of the surface gravity data used in carrying out the 
computations are described in Appendix A. 

3.1 Parameter Values Used for Computation 

In carrying out the initial geoid computations a value of W Q = 6263675.7 kgal 
m taken from Rap^ <1967) was used. The eilipscid parameters were given as 
1/f = 298.255 and a = 6378.155. The values of Ag s were made compatible with 
the chosen flattening by using AC 20 = .08 x 10" 6 and AC 40 = -.756 x 10“ 6 
when computing Ag s . The values of Ag e were made compatible by using a 
reference gravity formula of the form 

g = 978.0421(1 + 0.00530241 sin 2 q) - .0000059 sin 2 2cp) 



gure 1. Detailed Gravimetric Geoid of the United States 





Since it was known that the values of W 0 used was not completely correct 
it was anticipated that all the computed geoid undulation would be in error by a 
constant. 

In order to establish accurate absolute values of geoid undulation the 
geocentric station positions obtained using dynamic analyses were employed. 

3.2 Absolute Adjustment 

Table 2 presents a comparison of the geoid heights, computed in the present 
analysis, and geoid heights computed using the geocentric x, y, z positions ob- 
tained for 12 stations in North America in the SAO 1969 Standard Earth compu- 
tations (Gaposhkin and Lambeck, 1969). The SAO geoid heights are each given 
with reference to two ellipsoids, with semi-major axes of 6378.155, and 6378.137 
km, respectively; both with flattenings of 1/298.255. 

When the detailed gravimentric geoid heights at SAO station sites were 
compared to SAO geoid heights referenced to the 6378.155 km ellipsoid, a sys- 
tematic difference of 18 meters was noted. Rather than recomputing the gravi- 
metric geoid using an adjusted value for Wo, an equivalent adjustment was made 
to the value of a e characterising the reference ellipsoid to which the gravimetric 
geoid is assumed to be referenced. This consisted of subtracting 18 meters 
from the value originally used, which resulted in the value 6378.137 km reported 
above. 

As may be seen, the gravity geoid heights differ in a random manner by 10 
meters or less from the SAO station geoid heights referenced to an ellipsoid 


with a semi-major axis of 6378.137 km. 


rnmamm 


Table 2 


SAO - Gravity Geoid Comparisons 


( 1 ) 

Station 

No. 


(2) 

SAO Geoid* 
Height (155) 
(meters) 

(3) 

SAO Geoid 
Height (137) 
(meters) 

(4) 

Gravimetric 

Geoid 

(meters) 

-42 

-24 

-25 

-32 

-14 

-18 

-49 

-31 

-21 

-42 

-24 

-23 

-21 

- 3 

-11 

-42 

-24 

-24 

-52 

-34 

-29 

-35 

-17 

-21 

-39 

-21 

-26 

-42 

-24 

-24 

-46 

-28 

-19 

-37 

-19 

-29 


(3) - (4) 
(meters) 


Computed using data taken from Gaposhkin and Lambeck (1970) 
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Table 3 presents a similar comparison between the results of Marsh et al. 
(1970) and the gravimetric geoid. If one assumes Marsh's values are correct 
the geoid in Figure 1 should be considered as referenced to a geoid of semi- 
major axis 6378.124 kms. 

3.3 Comparative Evaluation 

In order to evaluate the accuracy of the final geoid, a number of compari- 
sons were made. The first comparison was made with the computations of 
Rapp (1967) of a combination geoid in a portion of the Western United States. 
Figure 2 presents a plot of these comparisons. In this compaxison, nine meters 
has been subtracted from the results of Rapp so that a comparison of relative 
shape is effected. No attempt at an absolute comparison was made. The agree- 
ment as to relative geoid shape is exceptional (differences exceeding 1,5 meters 
occur at only one point), consideriixg the entirely different sources of gravimetric 
data, and the differences in computational procedure. 

Another important source of comparative data is the detailed astrogeodetic 
geoid data. These data should have an rms accuracy on the order of 2 msters 
and should therefore provide an excellent basis for comparison throughout the 
United States. The most accurate astrogeodetic geoid information should be 
that of Rice (1970) which provide point astrogeodetic geoid values based largely 
on closed traverse loops with deflection measurements every 10 to 15 km. 
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Table 3 




Marsh - Gravity Geoid Comparisons 


(1) 

Station 

No. 

(2) 

Marsh Geoid* 
Height (155) 
(meters) 

(3) 

Marsh Geoid 
Height (124) 
(meters) 

(4) 

Gravimetric 

Geoid 

(meters) 

1021 

-55 

-24 

-25 

1022 

-50 

-19 

-18 

1030 

-53 

-22 

-28 

1034 

-53 

-22 

-18 

1042 

-57 

-26 

-21 

7037 

-60 

-29 

-23 

70^5 

-44 

-13 

-11 

7050 

-52 

-21 

•24 

7072 

-53 

-22 

-23 

7075 

-59 

-28 

-29 


Computed usinp data taken from Marsh et al. (1970) 


(5) 

(3) “(4) 
(meters) 
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Before any comparison can be made, the Rice Astrogeodetic geoid values 
must be transformed from the NAD to geoid values relative to a geocentric 
ellipsoid. This raises the questions of the transformation to use in converting 
NAD coordinates to geocentric coordinates. A number of transformations have 
been proposed, some involve simple transformations and others involve both 
transformations and rotations. Table 4 presents the differences between Rice’s 
astrogeodetic geoid and the gravimetric geoid after using each of four different 
sets of translation elements, and removing the mean differences. The impor- 
tant point to be noted is the overall degree of agreement between the astrogeo- 
detic and gravimetric geoid using the different translation coordinate sets. In 
all cases the rms differences are on the order of 2 meters. 

In addition to the point geoid heights. Rice (1967) has provided a detailed 
astrogeodetic geoid along the 35th parallel. Figure 3 graphically presents a 
comparison of this astrogeodetic geoid profile with the gravity geoid using the 
same procedure as employed with the point values and transformation set 3 of 
Table 4. Again the rms agreement is of the order of 2 meters, although at one 
point (95° longitude) the two differ by 4 meters. It is of interest to note, in 
Figure 2, that Rapp's result and the present gravimetric geoid agree well with 
one another at this point of maximum disagreement with the astrogeodetic 
value. 
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TABLE 


COMPARISON OF DETAILED GEOID vs RICE 1 S ASTROGEOOETIC GEOID UNDER 
VARYING ASSUMPTIONS FOR TRANSFORMING ASTROGEODETIC DATA 


Latitude 

Longi tude 

34° 

58' 

03V 0 

120° 

38' 

05V 5 

35 

00 

38.0 

119 

00 

48.0 

38 

47 

23.1 

121 

52 

15.6 

35 

02 

36.1 

106 

30 

24.1 

32 

13 

14.7 

106 

29 

41.6 

32 

00 

00.6 

103 

16 

07.2 

30 

59 

40.0 

098 

05 

50.5 

30 

36 

26.5 

091 

23 

18.1 

29 

38 

10.8 

091 

06 

49.3 

30 

59 

25.5 

089 

34 

29.5 

25 

30 

25.2 

080 

23 

17.5 

28 

29 

28.6 

080 

33 

35.6 

30 

36 

53.3 

081 

42 

14.8 

39 

28 

18.9 

076 

05 

15.2 

34 

59 

44.0 

076 

59 

11.7 

33 

28 

42.4 

091 

00 

08.5 

33 

34 

48.5 

092 

50 

07.2 

34 

56 

47.0 

093 

24 

18.3 

37 

38 

08.4 

094 

35 

46.8 

35 

03 

04.0 

097 

56 

52.6 

39 

13 

26.7 

098 

32 

30.5 

43 

37 

10.7 

096 

17 

52.3 

35 

06 

16.2 

103 

19 

55.0 

34 

56 

32.8 

096 

24 

55.3 

44 

43 

46.0 

105 

25 

50.7 

36 

47 

44.2 

103 

11 

48.5 

38 

50 

40.6 

102 

48 

46.8 

48 

06 

18.6 

102 

21 

09.7 

46 

44 

47.4 

102 

15 

13.4 

45 

12 

45.7 

102 

09 

14.1 

46 

21 

53.1 

108 

59 

07.3 

31 

03 

07.3 

102 

56 

05.8 
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Latitude 

Longitude 

| 

2 

3 

4 

41° 30' 41.9 

097° 37' 23.4 

1 

1 

1 

1 

30 48 49.8 

093 12 26.9 

-5 

-4 

-6 

-6 

47 50 28.9 

110 00 46.4 

1 

-1 

1 

1 


1* Corrected .fference between Rice's Astrogeodetic geoid and 
detailed gravimetric geoid. 

Using Marsh's translation values AX = - 25.1, AY = + 162.9, 
AZ s + 172.5 (Marsh, et.al., 1970) 

2* Corrected difference between Rice's Astrogeodetic geoid and 
detailed gravimetric geoid. 

Using Fischer's translation values AX = - 18 aY = + 145 
aZ * + 183 (Fischer, 1968) 

3* Corrected difference between Rice's Astrogeodetic geoid and 
detailed gravimetric geoid. 

Using SAO Standard Earth '66 translation values AX = - 30, 

AY = + 152, AZ = + 176. 

4* Corrected difference between Rice's Astrogeodetic geoid and 
detailed gravimetric geoid. 

Using SAO's translation values AX = - 25.8, aY = + 168.1 
AZ = + 167.0 (Lambeck, 1971) 







Figure 3. Detailed Gravimetric Geoid Vs Rice's Astrogeodetic Geoid Profile 











A final comparison was made using transformation sets 1 and 2 from 
Table 4 to transform geoid profile data taken from the map of Fischer (1967) to 
the geocentric system. Comparisons of the transformed astrogcodetic data and 
the gravimetric geoid along two East-West profiles located at latitude 40°N and 
45°N are presented in Figures 4 and 5. Again the relative agreement is within 
the ±2 meter range. 

4. CONCLUSIONS 

The conclusion which can be drawn, based upon comparisons made, is that 
the detailed gravimetric geoid presented here has precision of ±2 meters. 
Further study will be required to fix the absolute values of the geoid by choos- 
ing between the possibilities indicated in Tables 2 and 3 to provide an accuracy 
which is equivalent to the precision obtained. 

An important question which can be studied using the results presented 
here is the question of possible rotations of the North American Dataum. 

Figure 6 presents plots of the difference between the Astrogeodetic geoid 
heights of Rice which were transformed to geocentric-by a simple trans- 
lation and the detailed gravimetric geoid heights. The translation used was 
that of Marsh, i.e., AX = -25.1, Ay = +162.9, Az = +172.5. Examination of these 
plots indicate that there is no conclusive evidence of a rotation of the North 
American datum. If any datum rotation exists it is extremely small. Given 
the much smaller number of samples and the lesser accuracy available when 
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TRANSLATED ASTR0G30DETIC GEOID 
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-rttftf* 8 ’ 


deriving rotations using dynamically derived geocentric station positions, any 
rotations proposed on the basis of such data must be seriously questioned. 

5. ACKNOWLEDGMENTS 

The programming for this computation was carried out by Mrs. Brenda Hill 
while the actual computer runs were performed by Mr. Michael Nichols. 


22 


MM 




6. REFERENCES 

1. Bacon, D., Green, R., LeSchack, R., and Quintal, H., 1970, Geoid Represen- 
tation from Satellite Determined Coefficients, IBM Document No. TR-70-3. 

2. Fischer, I., 1967, Geoid Contours in North America from Astrogeodetic 
Deflections, 1927 North American Datum. Army Map Service. 

3. Fischer, I., 1968, A Modification of the Mercury Datum, Technical Report 
No. 67, U.S. Army Map Service, 53 pp. 

4. Gaposhkin, E. M., and Lambeek, K., 1970, 1969 Smithsonian Standard Earth 
(II), Special Report 315, Smithsonian Astrophysical Observatory, 93 pp. 

5. Lambeek, K., 1971, The Relation of Some Geodetic Datums to a Global Geo- 
centric Reference System, Bulletin Geodesique, No. 99, pp 37-54. 

6. Marsh, J. G., Douglas, B. C., and Martin, C. F., 1970, NASA STADAN, 
SPEOPT, and LASER Tracking Station Positions Derived from GEOS-1 and 
GEOS-2 Precision Reduced Optical and Laser Observations, in Proceedings 
of the GEOS-2 Program Review Meeting, Vol. 1, pp 187-214. 

7. Rapp, R. H:, 1967, Combination of Satellite and Terrestrial Data for a 
Detailed Geoid . Presented at the 48th Annual AGU Meeting April, 1967. 

8. Rice, D. A., 1970, Adjustment of Geoidal Sections in the United States, 

1927 North American Datum. Private Communication. 

9. Rice, D. A., 1962, A Geoidal Section in the United States , Bull. Geoid., 

Vol. 65 Pages 243-251. 

10. Vincent, S., Strange, W., Nichols, M., and Hill, B., 1971, Detailed Gravi- 
metric Geoid Computation, CSC Document No. 5023-05900-01 TR. 




I I'J ' JlitM 'fll 1 111 






APPENDIX A 

DERIVATION OF 1° X 1° TERRESTRIAL FREE- 
AIR ANOMALIES USED IN THE DETAILED GEOID COMPUTATION 

Except for portions of Canada, all of the terrestrial gravity data included 
herein were compiled from published sources. For many areas, particularly 
those outside continental margins, the data used were of a generalized nature, 
even though more detailed data are known to exist. These latter data were either 
unavailable or would require lengthy evaluation and analysis to obtain 1° x 1° 
mean anomaly values. Nontheless, all results are adjudged to be within the 
limits of accuracy required by the study. 

The following paragraphs include identification of all data sources, and 
description of computations used for conversion and averaging of various input 
data types. The areas of coverage of the different sources are shown in 
Figure A-l. 

A.l Areas 1, 2, and 3 - Continental Portions of U.S., Canada and Mexico 
Most of the coverage of these three areas is provided directly from 
previously published mean 1° x 1° free-air anomalies. Those of the continental 
U.S. (area 1) are from Strange and Woollard (1964); those of Mexico and Canada 
farea 2 and 3, respectively) are from Woollard (1968). A few 1° x 1° means 
were modified based on more recent information. 
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A. 2 Other Areas of Canada 


Values for the lighter hachured areas in Canada were derived from unpub- 
lished 2° x 1° mean free-air gravity anomalies, provided by Dr. D. Nagy, Gravity 
Division, Department of Energy, Mines aud Resources, Ottawa. Because of 
their generalized nature, these values were applied only where more detailed 
data were lacking, or were suspect. 

In most places, the 2° x 1° values were applied directly without modification. 
However, where one of the 1° x 1° squares within the 2° x 1° areas was previously 
defined by a more accurate data source, the value of the blank 1° x 1° was de- 
fined such that the average for it and the previous value equalled that of the 
2° x 1°. 


A.3 Area 4 - Washington - Oregon Coast. 

The data sources for this area were free-air gravity contour maps pub- 
lished by Dehlinger (1969). Unless the setting of a 1° x 1° square was par- 
ticularly straight forward, it was subdivided into segments of approximately 
equal gravity. The mean value for the 1° x 1° square then, was computed as the 
average of the mean values of the individual segments, weighted proportionately 
relative to the segment area. 

A. 4 Area 5 - Central California Coast 

Free-air anomalies for 1° x 1° squares of area 5 were derived from the 
Bouguer and bathymetric map of Lattimore, Bush and Bush (1968). Each 





1° x 1° square of their map was divided into nine equal areas, from which the 
Bouguer gravity and water depth of the midpoint of each of the nine segments 
were approximated from the source maps. A free-air value at each midpoint 
was then computed by the following equation from Nettleton (1940, p. 54) 

F.A. = B.A. - (.04185) ( o-)(h) milligal/meter 


where 

F.A. = Free-air anomaly 
B.A. = Bouguer anomaly 

h = water depth below sea level, meters 
cr = specific gravity 
= 2.84 - 1.03 3 
= 1.81 gm/cm 3 

The latter density values are those applied by Lattimore et al. to derive the 
Bouguer anomalies. 

The mean free-air anomaly of each 1° x 1° square is the average of its 
nine component free-air values. 

A.5 Area 6 - Baja California Coast 

The free-air values of 1° x 1° squares within area 6 were derived by 
plotting and contouring all free-air measurements within the area as reported 
by Worzel (1965). Mean values for 1° x 1° squares were approximated from 
the 10 milligal contours by methods described for area 4. 
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A. 6 Area 7 - Gulf of Mexi co 

Values for 1° x 1° squares of this area were derived from the free-air 
gravity map of Dehlinger and Jones (1965), in the same manner as that applied 
for the area 4 data. 

A. 7 Area 8 - Florida 

The free-air gravities of the ten 1° x 1° squares of Florida were estimated 
directly from Woollard's Bouguer Gravity Map of the United States (1964), by 
methods described for area 4. 

A. 8 Area 9 

The free-air values of area 9 were derived from three difference sets of 
data: 

(a) Free-air gravity measurements of Cuba, mostly located ir the northern 
and eastern portions, as published by Dickerson (1940), and Shubert 
(1957). 1° x 1° squares for these portions of Cuba were visually esti- 
mated from 10 milligal contours of these data, as supplemented by off 
shore data, mentioned below. 

(b) A Bouguer map of Southern Andros Island and shallow adjacent marine 
areas, computed by Richards and Malone (1949). Because of the low 
elevations on Andros Island and the shallowness of the sea over a broad 
area south of the island, test computations showed that differences be- 
tween Bouguer and free-air gravity is less than a milligal. Therefore, 
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formal conversion was unwarranted, and 1° x 1° values for this area 
were visually estimated directly from the Bouguer contour map. 

(c) Free-air values at sea, provided by Worzel (1965). These data were 
relatively sparce considering the topographic variations of the sea 
floor and island chains, and hence, could not be dependably contoured. 
Rather, where a number of readings were listed within the same gen- 
eral area, and apparently shared similar conditions including water 
dpeth, their arithmetic average was assumed to represent the general 
area as well as adjacent areas of similar characteristics where data 
was lacking. Values derived from areas 7 and 10 were similarly 
extrapolated into data- poor segments of the Andros- Cuba area. 

A. 9 Area 10 - Continental Shelf of Southeastern U.S. 

1° x 1° mean values were interpolated directly from the free-air gravity 
map of the eastern U.S. continental shelf, contoured at 20 milligal intervals by 
Emory, et al (1970). The mean value for each 1° x 1° square was visually 
estimated on the basis of proportional areas as for area 4. In areas with mini- 
mal gravitational relief, the method was aided by adding supplemental 10 
milligal contours. 

A. 10 Area 11 - Geo rgia Coast 

The mean values of the three 1° x 1° squares of this area were estimated 
directly from the free-air 5 milligal contour map published by U.S. Coast and 
Geodetic Survey (1968). 
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A. 11 Area 12 - Contine ntal Shelf of Northeastern U.S. 

Data for tills area were free-air values provided by Worzel (1965), for 
various measurements prior to 1959, and those derived from Chain Cruises 70 
and 73, published by Bowen and Aldrich (1969a, 1969b). Data from all three 
sources were plotted together and cuntoured at 10 milligal intervals. Values 
for each 1° x 1° within the area we: o then estimated by methods applied for 
area 4. 

A. 12 Area 13 - Gulf of Maine 

Five of the squares vuthin this area were estimated directly from the 5 
milligal contours of the free-air gravity map prepared by Yellin (1968). The 
value for the sixth square in the northeast corner of the area, where contours 
were generally lacking, was the arithmetic average of the 1° x 1° free-air 
values within the square, as provided by Yellin. 

A. 13 Area 14 - The Gulf of St. Lawrence 

The data sources for area 14 were bathymetric depths and simple Bouguer 
anomalies provided by Goodacre, Brule, and Cooper (1969). The Bouguer values 
were converted into free-air values by applying the same equation as for area 
5, except that rock density was 2.67. These values were supplemented by 
Bouguer gravity along the major coast lines, shown by the 50 milligal controur 
lines of the gravity map of Goodacre, et al. 

Results were plotted and contoured at 10 milligal intervals, allowing visual 
estimates of mean free-air values for each of the 1° x 1° squares within the area. 
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A. 14 Area 15 - Northern Lake Superior 

Free-air anomalies for area 15 were computed from the Bouguer gravity 
map of Lake Superior by Weber and Goodacre (1966), in conjunction with the 
bottom topography furnished by Wold and Ostenso (1966). Each 1° x 1° square 
of the area on each map was divided into 16 equal portions, and the water depth 
and Bouguer gravity at each line intersection were tabulated. Water depths were 
converted into elevations assuming a lake elevation of 166 feet above sea level. 
Then free-air values were computed for each intersection by the equation, 

F.A. = B.A. + B.C. 

where 

F.A. = Free-air anomaly 

B.A. = Bouguer anomaly 

B.C. = Bouguer correction = a + b - c 

a = correction for water above sea level 
= feet of water above sea level x water density x m 
b = correction for rock above sea level 
= feet of rock above sea level x rock density x m 
c = correction for water below sea level 
= (feet of water below sea level) (rock density - water density) x (m) 
m = constant, .01276, after Nettleton (1940, p. 54) 

Rock and water densities are 2.67 and 1.00 respectively, as assumbed by Weber 
and Goodacre. 
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